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We use QCD sum rules to study the possible existence of QQ — ud mesons, assumed to be a 
state with J p = 1 + . For definiteness, we work with a current with an axial heavy diquark and a 
scalar light antidiquark, at leading order in a s . We consider the contributions of condensates up 
to dimension eight. For the &-quark, we predict Mr bb = (10.2 ± 0.3) GeV, which is below the BB* 
>— ^ ' threshold. For the c-quark, we predict Mt cc = (4.0 ± 0.2) GeV, in agreement with quark model 

predictions. 
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The general idea of possible stable heavy tetraquarks has been first suggested by Jaffe pj . The case 
of a tetraquark QQud with quantum numbers I — 0, J — 1 and P = +1 which, following ref. 0, we call 
Tqq, is especially interesting. As already noted previously 0, Q, the Tf>b and T cc states cannot decay 
strongly or electromagnetically into two B or two D mesons in the S wave due to angular momentum 
conservation nor in P wave due to parity conservation. If their masses are below the BB* and DD* 
f-- thresholds, these decays are also forbidden. Moreover, in the large tuq limit, the light degrees of freedom 

cannot resolve the closely bound QQ system. This results in bound states similar to the Aq states, with 
QQ playing the role of the heavy antiquark 0]. Therefore, the stability of Aq implies that QQud is also 
safe from decaying through QQud — > QQq + qud . As a result, Tqq is stable with respect to strong 
interactions and must decay weakly. 

There are some predictions for the masses of the Tqq states. In ref. [5[ the authors use a color-magnetic 
interaction, with flavor symmetry breaking corrections, to study heavy tetraquarks. They assume that 
the Belle resonance, X(3872), is a cqcq tetraquark, and use its mass as input to determine the mass of 
other tetraquark states. They get M Tcc = 3966 MeV and M Tbb = 10372 MeV. In ref. 2], the authors 
use one-gluon exchange potentials and two different spatial configurations to study the mesons T cc and 
T bb . They get M Trc = 3876 - 3905 MeV and M Thb = 10519 - 10651 MeV. There are also calculations 
using expansion in the harmonic oscillator basis [SL and variational method 

In this work we use QCD sum rules (QCDSR) [H, H, 03, to study the two-point functions of the state 
Tqq. There are several reasons, why it is interesting to investigate this channel. First of all, having 
two heavy quarks, it is an explicit exotic state. The experimental observation would already prove the 
existence of the tetraquark state without any theoretical extrapolation. Moreover, from a technical point 
of view, this means that there are no contributions from the disconnected diagrams, which are technically 
very difficulty to estimate in QCD sum rules or in lattice gauge theory calculation. 

Inprevious calculations, the QCDSR approach was used to study the light scalar mesons [III E3, EH 
EH, EH the D^"j(2317) meson 0, [TtJ and the X(3872) meson [T8[, considered as four-quark states and a 
good agreement with the experimental masses was obtained. However, the tests were not decisive as the 
usual quark-antiquark assignments also provide predictions consistent with data [ToL [l2l [l^ . |20| . 

Considering Tqq as an axial diquark-antidiquark state, a possible current describing such state is given 
by: 

U = i[Qlc~tMu a ^cdJ], (i) 

where a, b are color indices, C is the charge conjugation matrix and Q denotes the heavy quark. 

In general, one should consider all possible combinations of different 1 + four-quark operators, as was 
done in [2l[ for the ++ light mesons. However, the current in Eq.|T|) well represents the most attractive 
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configuration expected with two heavy quarks. This is so because the most attractiv e lig ht antidiquark 
is expected to be the in the color triplet, flavor anti-symmetric and spin channel |22l l23l. [2~j |. This 
is also expected quite naturally from the color magnetic interaction, which can be phcnomcnologically 
parameterized as, 

C 

Vij — Aj • XjOi ■ <7j. (2) 

J m l m,j 

Here, to, A, a are the mass, color and spin of the constituent quark Eq.© favors the anti-diquark to 
be in the color triplet and spin channel. The flavor anti-symmetric condition then follows from requiring 
anti-symmetric wave function of the anti-diquark. Similarly, since the anti-diquark is in the color triplet 
state the remaining QQ should be in the color anti-triplet spin 1 state. Although the spin 1 configuration 
is repulsive, its strength is much smaller than that for the light diquark due to the heavy charm quark 
mass. Therefore a constituent quark picture for Tqq would be a light anti-diquark in color triplet, flavor 
anti-symmetric and spin (e a b c [ub'fs C3^\) combined with a heavy diquark of spin 1 (e ae f[Q^Cj^Qf]). 
The simplest choice for the current to have a non zero overla p w ith such a Tqq configuration is given 
in Eq. (Q]). While a similar configuration T ss is also possible (25|, we believe that the repulsion in the 
strange diquark with spin 1 will be larger and hence energetically less favorable. As discussed above, 
since the quantum number is 1 + , the decay into DD or BB would be forbidden and the allowed decay 
into DD* or BB* would have a smaller phase space, and the tetraquark state might have a small width, 
or may even be bound. 

The QCDSR is constructed from the two-point correlation function 

U^(q) = ijd±x e^(0\T[j,(x)jt(0)}\0) = -I^g 2 )^ - Sip.) + H (q 2 )^. (3) 

Since the axial vector current is not conserved, the two functions, ITi and ITo, appearing in Eq. ([3]) are 
independent and have respectively the quantum numbers of the spin 1 and mesons. 

The calculation of the phenomenological side proceeds by inserting intermediate states for the meson 
Tqq. Parametrizing the coupling of the axial vector meson 1 + , to the current, j M , in Eq. |T]) in terms of 
the meson decay constant fx and the meson mass Mp as: 

(0\MTqq) = V2f T M^ , (4) 
the phenomenological side of Eq. © can be written as 



<"'<' 2 ) = f^ (-■'- + if) 



where the Lorentz structure gets contributions only from the 1 + state. The dots denote higher 
axial-vector resonance contributions that will be parametrized, as usual, through the introduction of a 
continuum threshold parameter so (26j . 

On the OPE side, we work at leading order in a s and consider the contributions of condensates up to 
dimension eight. To keep the charm quark mass finite, we use the momentum-space expression for the 
charm quark propagator. We follow ref. [13] and calculate the light quark part of the correlation function 
in the coordinate-space, which is then Fourier transformed to the momentum space in D dimensions. The 
resulting light-quark part is combined with the charm quark part before it is dimensionally regularized 
at D = 4. 

The correlation function, IIi, in the OPE side can be written as a dispersion relation: 

n? pf V) = r ds I^L t (6) 

J 4 m 2 Q s-r 

where the spectral density is given by the imaginary part of the correlation function: 7r/?(s) = 
Im[nf >ps (s)]. After making a Borel transform of both sides, and transferring the continuum contri- 
bution to the OPE side, the sum rule for the axial vector meson Tqq up to dimension-eight condensates 
can be written as: 

2/2 M |e- M -/ M2 = f ° ds e- s ' M * p(s) +nr x{ciq) (M 2 ) , (7) 
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where 



with 



p(s) = pP ert (s) + {s) + p V> (s) + p mix (s) + p™' (s) + p 



(qq) 2 



rnix(qq) . 
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(9) 



where the integration limits are given by a m i n = (1 — ^/l — 4jtoq/s)/2, a max = (1 + «/l — 4toq/s)/2 

and /3 m j n = ar7ig/(sa — TOg). The contribution of dimension-six condensates (g 3 G 3 ) is neglected, since it 
is assumed to be suppressed by the loop factor 1/167T 2 . We have included, for completeness, a part of the 
dimcnsion-8 condensate contributions. We should note that a complete evaluation of these contributions 
require more involved analysis including a non-trivial choice of the factorization assumption basis [28| 
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In order to extract the mass Mt without worrying about the value of the decay constant fx, we take 
the derivative of Eq. ([7]) with respect to 1/M 2 , divide the result by Eq. ([7]) and obtain: 



Aff. 



JZ 2 ds e-°/M* p(s) ■ 



(11) 



This quantity has the advantage to be less sensitive to the perturbative radiative corrections than the 
individual moments. Therefore, we expect that our results obtained to leading order in a s will be quite 
accurate. 

In the numerical analysis of the sum rules, the values used for the quark masses and condensates are (see 
e.g. 0,111): m c (m c ) = (1.23 ±0.05) GeV, m b {m b ) = (4.24 ±0.06) GeV, (qq) = -(0.23 ± 0.03) 3 GeV 3 , 
(qga.Gq) = m%{qq) with m 2 = 0.8 GeV 2 , (g 2 G 2 ) = 0.88 GeV 4 . 

We start with the double charmed meson T cc . We evaluate the sum rules in the range 2.0 < M 2 < 
4 GeV 2 for s in the range: 4.6 < ^ < 5.0 GeV. 

Comparing the relative contribution of each term in Eqs. © to (TIT)]) , to the right hand side of Eq. {7J 
we obtain a quite good OPE convergence (the perturbative contribution is at least 50% of the total) 
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FIG. 1: The relative OPE convergence in the region 2.0 < M 2 < 4.0 GeV 2 for ^/sE = 4.8 GeV. We start with the 
pcrturbative contribution divided by the total (long-dashed line) and each subsequent line represents the addition 
of one extra condensate dimension in the expansion: + {g 2 G 2 } (dot-dashed line), + {qq) 2 (dotted-line), +m 2 (qq} 2 
(solid line). 



for M 2 > 2.5 GeV 2 , as can be seen in Fig. Q] This analysis allows us to determine the lower limit 
constraint for M 2 in the sum rules window. This figure also shows that, although there is a change of 
sign between dimension-six and dimension-eight condensates contributions, the contribution of the latter 
is very small, where, we have assumed, in Fig. [T] to Fig. 2J the validity of the vacuum saturation for 
these condensates. The relatively small contribution of the dimension-eight condensates may justify the 
validity of our approximation, unlike in the case of the 5-quark current correlator, as noticed in [301 ] . 
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FIG. 2: The solid line shows the relative pole contribution (the pole contribution divided by the total, pole plus 
continuum, contribution) and the dashed line shows the relative continuum contribution for ^/sq = 4.8 GeV. 



We get an upper limit constraint for M 2 by imposing the rigorous constraint that the QCD continuum 
contribution should be smaller than the pole contribution. The maximum value of M 2 for which this 
constraint is satisfied depends on the value of so- The comparison between pole and continuum contri- 
butions for ^/Jq — 4.8 GeV is shown in Fig. O The same analysis for the other values of the continuum 
threshold gives M 2 < 3.1 GeV 2 for = 4.6 GeV and M 2 < 3.6 GeV 2 for = 5.0 GeV. 

In Fig. [3l we show the T cc meson mass obtained from Eq. (fTTj) , in the relevant sum rules window, with 
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FIG. 3: The T cc meson mass as a function of the sum rule parameter (M 2 ) for different values of the continuum 
threshold: y/sio = 4.6 GeV (dotted line) and y/so = 5.0 GeV (solid line). The bars indicate the region allowed 
for the sum rules: the lower limit (cut below 2.5 GeV 2 ) is given by OPE convergence requirement and the upper 
limit by the dominance of the QCD pole contribution. 



the upper and lower validity limits indicated. From Fig. [3] we see that the results are reasonably stable 
as a function of M 2 . In our numerical analysis, we shall then consider the range of M 2 values from 2.5 
GeV 2 until the one allowed by the sum rule window criteria as can be deduced from Fig. [3] for each value 
of s . 

We found that our results are not very sensitive to the value of the charm quark mass, neither to 
the value of the condensates. The most important source of uncertainty is the value of the continuum 
threshod and the Borel interval. Using the QCD parameters given above, the QCDSR predictions for the 
T cc mesons mass is: 

M Tcc = (4.0 ± 0.2) GeV, (12) 

in a very good agreement with the predictions in refs. @ and @. 

One can also evaluate the decay constant, defined in Eq. (HJ), to leading order in a s : 

f Tcr = (5.95 ± 0.65) x 10" 5 GeV , (13) 

which can be more affected by radiative corrections than Mj ce . 

In the case of the double-beauty meson Tbb, using consistently the perturbative MS'-mass m^irrn,) = 
(4.24±0.6) GeV, and the continuum threshold in the range 11.3 < y/so < 11.7 GeV, we find a good OPE 
convergence for M 2 > 7.5 GeV 2 . We also find that the pole contribution is bigger than the continuum 
contribution for M 2 < 9.6 GeV 2 for y/aE < 11.3 GeV, and for M 2 < 11.2 GeV 2 for ^ < 11.7 GeV. 

In Fig. 0] we show the Tbb meson mass obtained from Eq. (jTTjl , in the relevant sum rules window, with 
the upper and lower validity limits indicated. From Fig. [4] we see that the results are very stable as a 
function of M in the allowed region. Taking into account the variation of M 2 and varying s and m,b in 
the regions indicated above, we arrive at the prediction: 

M Tbb = (10.2 ± 0.3) GeV , (14) 

also in a very good agreement with the results in refs. 0, @ and 0. For completeness, we predict the 
corresponding value of the decay constant to leading order in a s : 

f Tbb = (10.4 ± 2.8) x TT 6 GeV . (15) 

We have presented a QCDSR analysis of the two-point functions of the double heavy-quark axial meson, 
Tqq, considered as a four quark state. We find that the sum rules results for the masses of T cc and Tbb 
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FIG. 4: The Tj,& meson mass as a function of the sum rule parameter (M 2 ) for different values of the continuum 
threshold: ^/so = 11.3 GeV (dashed line), y^So = 11.7 GeV (solid line). The bars delimit the region allowed for 
the sum rules. 

are compatible with the results in refs. [2| and [5[. An improvement of this result needs an accurate 
determination of running masses m c and of the MS-scheme and the inclusion of radiative corrections. 

Our results show that while the T cc mass is bigger than the D*D threshold at about 3.875 GeV, the 
Tbb mass is appreciably below the B*B threshold at about 10.6 GeV. Therefore, our results indicate that 
the Tbb meson should be stable with respect to strong interactions and must decay weakly. Our result 
also confirms the naive expectation that the exotic states with heavy quarks tend to be more stable than 
the corresponding light states(3lj|. 

We present in Eqs. (JT3J), and (fT5|) predictions for the decay constants of the T cc and Tbb- 
Different choices of the four-quark operators have been systematically presented for the ++ light 
mesons in [2l|. Though some combinations can provide a faster convergence of the OPE, we do not 
expect that the choice of the operators will affect much our results, where, in our analysis, the OPE has 
a good convergence. 
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